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The  main  objectives  of  the  original  proposal  were  to 
calculate  spontaneous  and  other  polarization  properties  of 
polar  and  piezoelectric  crystals.  The  purpose  of  the  work 
was  to  aid  in  selecting  the  best  pyroelectric  materials  to  be 
used  in  pyroelectric  image  tubes.  The  crystals  for  which 
the  main  work  effort  was  directed  were:  NaNtf^^  LiNbO^ 


and  KNbO^.  Calculations  on  the  first  three  crystals  have 
led  to  five  publications  and  calculations  on  the  last  are 
being  compiled  for  publication.  Publications  on  other  aspects 
of  the  grant  are  listed  in  the  List  of  Publications  which 
includes  the  five  publications  mentioned  above. 

The  computer  piragram  used  to  make  the  calculations 
was  written  by  R.  N.  ^wema  formerly  at  the  Aerospace  Research 
Laboratory,  WPAFB,  Ohio  reusing  charge  densities  obtained  from 
band  calculations  with  this  program,  the  following  piezoelectric 
coefficients  <32i»  ^22'  ^23  sodium  nitrite  have  been 

calculated.  This  was  accomplished  by  calculating  the  changes 
in  the  electric  polarization  of  sodium  nitride  for  three 
different  strains.  The  piezoelectric  coefficients  are  linear 
relations  relating  induced  polarizations  to  applied  stresses. 
Experimental  elastic  coefficients  were  used  to  obtain  stresses 
from  the  above  strains.  There  are  some  vmcertainties  in  the 
experimental  elastic  coefficients  ^23* 


The  uncertainties 
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in  the  coefficients  ^23  rise  to  uncertainties  in 

the  piezoelectric  coefficients  <322^  and  <323«  Our  calculated  value 
for  the  piezoelectric  coefficient  compares  well  with  the 
experimental  value  while  the  calculated  values  of  d2jL  uud  d23 
agree  only  fairly  well  with  experimental  values.  The  details 
of  the  calculations  and  the  results  of  the  calculations  are 
presented  in  the  article  of  reference  (2) . 

Nonlinear  strain  -  polarization  calculations  were  made 
on  crystalline  NaN02.  These  calculations  showed  a  dependence 
of  polarization  of  strain  that  had  a  discontinuous  slope.  A 
short  article  on  a  possible  second  order  phase  transition  in 
crystalline  NaN02  due  to  strain  is  to  be  pviblished  in  Physics 
Letters. 

The  plain-wave-Gaussian  mixed  basis  method  has  been 
employed  to  calculate  crystalline  energy  bands  for  cubic 
boron  nitride (zinc-blende  structure).  These  calculations  agree 
favorably  with  experiment.  Polarization  calculations  for 


sheared  piezoelectric  cubic  boron  nitride  have  been  undertaken. 
With  a  shear  strain  £4  =  .001  along  the  boron-nitrogen  bond 
which  is  directed  along  the  body  diagonal  in  the  unit  cell,  the 
polarization  is  changed  about  ”9.56  x  10  e.s.u./cm  .  In 
order  to  calculate  the  corresponding  piezoelectric  constant 
the  stiffness  constants  for  boron  nitride  are  needed.  So  far 
these  have  not  been  found.  However,  they  can  be  estimated  to 
determine  piezoelectric  coefficients. 
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The  spontaneous  polarization  of  LiNbO^  has  been  calculated 
and  the  value  differs  by  20  percent  from  the  experimental 

•  ff,  f  ^  '  ■  '■  ^  r 

value.  Because  contains  so  many  electrons  in  a  unit 

cell  we  were  not  able  to  do  a  thorough  calculation  on  LiNbO^ 
as  was  done  on  NaNd^. 

The  article  entitled  "Off-Diagonal  Disorder  in  Dilute 
Alloys"  has  been  published  in  Physical  Review  B.  The  research 
on  the  work  was  sponsored  in  part  by  the  grant.  The  content 
of  the  paper  is  to  resolve  the  difference  between  the  results 
of  Fukuyama,  Krakauer,  Schwartz  and  Hwang  and  Sher  on  the 
electronic  density  of  states  in  dilute  metal  alloys  with  off- 
diagonal  disorder  retained  in  the  Hamiltonia.  The  work  is 
accon^lished  by  introducing  an  alternative  decomposition  of 
the  alloy  Hamiltonia  in  the  Hwang  and  Sher  fonnalism. 

The  Operator  Recursion  Method  was  previously  applied 
to  the  motion  of  a  single  particle  moving  in  a  1-dimensional 
periodic  lattice.  Work  on  the  method  in  this  final  year 
of  the  grant  involved  its  extension  to  3-spa tial  dimensions 
and  to  the  meiny  body  problem.  The  previous  cumbersome 
procedure  of  the  self-consistent  removal  of  momentum 
operators  from  the  recursion  equations  has  been  replaced  by 
a  generating  function  technique  '  .  The  recursion  relations 
are  now  developed  around  a  matrix  element  of  products  of 
a  periodic  functions  of  position  with  powers  of  the  momentum 
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operator.  The  recursion  relation  has  been  transformed 
into  a  linear  partial  differential  equation  for  the 
generating  function  having  a  number  of  independent  variables 
equal  to  the  number  of  degrees  of  freedom,  plus  time.  The 
charge  density  for  the  system  is  determined  by  the  initial 
conditions  of  the  generating  function,  while  the  band 
structure  follows  from  seperating  out  the  time  variable 
and  solving  the  associated  eigenvalue  problems.  The  problem 
of  the  separation  of  variables  is  mitigated  by  group 
theoretical  results  following  from  the  symmetries  of  the 
system. 

The  generating  function  procedure  is  also  employed 
ip  approaching  the  many  body  problem.  Here  the  recursion 
relation  is  developed  for  matrix  elements  of  normal  products 
of  raising  and  lowering  operators.  For  a  1-dimensional 
system  involving  only  coupling  between  the  "charge”  density 
with  an  external  potential  the  generating  function  involves 
3  time  and  2  additional  independent  variables.  Again  the 
band  structure  and  charge  density  follow  from  the 
separation  of  variables  and  the  linear  PDE  initial  condi¬ 
tions.  More  complicated  systems  involving  self  coupling 
of  the  field  lead  to  an  increase  in  the  number  of  independent 
variables,  the  number  being  a  function  of  the  number  of 
degrees  of  freedom.  Relaxation  of  the  system  is  described 
by  solutions  of  the  PDE  which  are  associated  with  non-self 
adjoint  operators  acting  on  the  generating  function. 
Investigations  are  continuing  on  the  classification  and 
application  of  these  equations. 
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The  interaction  of  a  ferro  or  pyroelectric  system  with 
an  external  field  in  the  dipole  approximation  was  previously 
described  (in  our  Proposal  1977) .  The  polarization  relax¬ 
ation  fvinction  describes  how  such  a  system  responds  to  an 
external  perturbation.  It  is  also  essential  for  calculating 
other  time-dependent  properties  of  ferro  and  pyroelectric 
materials.  Taking  advantage  of  the  two  recent  methods 
developed  by  us,  we  have  calculated  the  relaxation  function 
for  a  van  der  Waals  charge  system.  A  charge  density  system, 
e.g.,  a  ferroeletric  material,  where  all  charges  interact 
with  all  others  with  equal  strength  is  called  a  van  der 
Waals  charge  system.  Although  this  model  system  greatly 
simplies  many  complexities  of  real  systems,  it  is  now  known 
that  this  model  treats  static  and  dynamic  properties 
sensibly  as  long  as  the  systems  are  not  in  the  critical 
region , 

The  relaxation  function  is  defined  by 
R(t)  =  Z”^  Tr{e“^”  P(t)  P} 

where  Z  is  the  partition  function,  P  is  the  charge  density 
operator.  Using  the  van  der  Waals  charge  model  for  H,  we 
have  deduced  analytic  expressions  for  R(t)  in  different 
spatial  symmetry  regimes.  (Previously  we  had  obtained  only 
leading  factors.  We  now  have  complete  solutions.)  The 
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spatial  symmetry  enters  into  the  problem  through  spatial 
orientation  of  the  interaction  given  in  H.  If  the  inter¬ 
action  is  spatially  random,  one  has  a  Heisenberg-like 
rotational  symmetry.  If  the  interaction  is  dominantly 
planar,  one  has  an  XY-like  rotator  symmetry.  If  the 
interaction  is  linear,  one  has  an  Ising-like  asymmetry. 

We  have  now  obtained  exact  solutions  for  a  number  of  pure 
symmetries  and  mixed  symmetries. 

Here  we  shall  summarize  our  results  wbich  will  appear 
in  Physical  Review  1979.  Our  general  results  fall  into  two 
physical  regimes:  above  and  below  the  transition  temperature 
If  temperature  is  above  the  transition  temperature,  the  relax 
ation  function  is  exactly  Gaussian  for  both  the  linear 
and  planar  spatial  symmetries, 

R(t)  =  R(0)  exp- (At) ^ 

where 

R(0)  =  <P^>  . 

The  relaxation  function  for  the  rotational  symmetric  case 
is  uninteresting  as  it  vanishes  with,  the  vanishing  of  the 
perturbation.  This  arises  from  random  orientation  of 
interaction  which  does  not  discriminate  the  manner  in  which 
it  is  perturbed.  For  the  other  pure  symmeric  cases,  the 
widths  of  the  Gaussian  A  ^  are  of  importance.  For  the 
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planar  spatial  symmetry,  the  width  is  constant,  not  depending 
on  temperature.  Hence,  the  relaxation  function  "dies" 
very  smoothly  in  time  unaffected  by  fluctuations  in  charge 
densities  brought  about  by  thermal  energies.  For  the  pure 
linear  spatial  symmetry,  the  width  is  temperature-dependent. 
Hence,  the  relaxation  function  for  this  case  does  not  behave 
monotonically.  As  the  temperature  is  lowered  toward  the 
critical  temperature  for  example,  there  is  narrowing  of 
the  width  resulting  in  critical  slow  down  of  all  relaxation 
processes . 

If  temperature  falls  below  the  transition  temperature, 

the  shape  of  the  relaxation  function  now  becomes  dependent 

on  the  nature  of  spatial  symmetries  directly.  For  the 

pure  planar  spatial  symmetry  the  relaxation  function  shape 

remains  Gaussian  whose  width,  however,  is  no  longer  temper- 

—1/2 

ature-independent.  It  has  a  T  '  -dependence  in  such  a 
way  that  it  grows  significantly  with  temperature.  For  the 
pure  linear  spatial  symmetry,  the  relaxation  function  is 
oscillatory, 

R(t)  =  R(0)cos  wt  +  Q(0)  sin  wt, 

where  Q  is  to  a  leading  order  a  time  independent  constant 
depending  on  temperature,  the  frequency  is  a  precessional 
frequency  depending  on  the  degree  of  internal  ordering. 

It  vanishes,  for  example,  at  the  critical  temperature 
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suffering  critical  slow  down. 

The  advantage  of  using  the  van  der  Waals  model  for 

complex  systems  becomes  apparent  from  these  solutions,  in 

particular  from  the  correspondence  we  have  established  for 

one  dimensional  problems.  In  the  high  temperature  limit 

(i.e.,  T-*^) ,  numerical  and  analytical  solutions  for  1- 

dimensional  nearest-neighbor  interaction  model  indicate 

2  2 

that  R(t)'“exp(-J  t  )>  where  J  is  the  strength  of  the  nearest 
neighbor  charge  coupling.  This  result  is  directly  obtain¬ 
able  from  our  van  der  Waals  solutions  if  one  properly 
adjusts  the  range  of  interaction  in  the  H  used  for  the  van 
der  Waals  model.  Thus,  even  for  highly  complex  ferro 
and  pyroelectric  materials,  their  time-dependent  behavior 
is  qualitatively  correctly  described  by  the  van  der  Waals 
model.  (Estciblishing  this  result  was  the  main  purpose  of 
this  phase  of  our  proposed  work) .  For  quantititive  accuracy 
one  must  forgo  the  analytic  approach  and  utilize  high-speed 
computers.  This  may  be  done  by  adding  to  the  van  der  Waals 
mode  some  additional  features  of  interaction  but  preserving 
the  spatial  symmetries.  These  results  and  their  applica¬ 
tions  to  other  dynamical  problems  have  been  reported  at 
a  recent  international  physics  conference  held  in  Grenoble 
France.  The  details  of  this  contribution  will  appear  in 
the  Journal  de  Physique  (paris) . 
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